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TRANSIENT  RESPONSE  OF  LAMINATED,  BIMODULAR-MATERIAL, 
COMPOSITE  RECTANGULAR  PLATES 


J.  N.  Reddy f 

Virginia  Polytechnic  Institute  and  State  University 
Blacksburg,  VA  24061 


The  paper  presents  a  finite-element  analysis  of  the  transient 
behavior  of  fiber-reinforced,  single-layer  and  two-layer  cross-ply 
rectangular  plates  of  bimodular  materials  (i.e.,  materials  whose  linear 
elastic  properties  are  different  depending  on  whether  the  fiber- 
direction  strains  are  tensile  or  compressive).  To  validate  the  finite 
element  results,  a  closed-form  solution  is  also  presented  for  a 
rectangular  plate  with  all  edges  simply-supported  without  in-plane 
restraint  (along  the  edge)  and  tangential  rotation  and  subjected  to 
suddenly  applied  sinusoidally  distributed  normal  pressure.  The  time 
behavior  of  the  transverse  loading  is  arbitrary  (e.g.,  step  loading, 
impulse  loading,  etc.).  Numerical  results  for  transverse  deflection  and 
locations  of  the  'neutral  surface1  as  functions  of  time  are  presented 
for  two  bimodular  materials.  The  finite  element  results  agree  very 
closely  with  the  closed-form  solutions. 


Introduction 

In  1941  S.  Timoshenko  [1]  considered  the  flexural  stresses  in 
materials  which  have  different  moduli  in  tension  and  in  compression. 
Examples  of  such  materials,  called  bimodular  materials,  are  provided  by 
cord-rubber  composites,  certain  biological  tissues,  and  paperboard, 
among  others.  Analysis  of  bimodular  (or  bimodulus)  materials  in  two 
dimensions  began  with  the  work  of  Ambartsumyan  [2]  in  1965.  Following 
this  work,  there  appeared  numerous  static,  two-dimensional  analyses  of 
bimodular  materials.  Among  these,  the  works  of  Kamiya  [3-5],  Jones  and 
Morgan  [6],  and  Bert,  Reddy  and  their  colleagues  [7-11]  should  be 
particularly  noted  in  the  context  of  bending  of  bimodular  plates.  For 
additional  references  to  the  subject,  the  reader  is  referred  to  the 
bibliography  in  [7-11].  To  the  best  of  the  present  author's  knowledge, 


^Professor,  Department  of  Engineering  Science  and  Mechanics. 

i 


2 


no  previous  transient  analyses  of  plates  laminated  of  bimodular 
composite  materials  are  available  in  the  open  literature.  The  present 
work  is  believed  to  be  the  first  one  to  consider  the  transient  analysis 
of  layered  composite  plates  of  bimodular  materials. 

The  present  study  employs  the  finite  element  model  developed  in  [9] 
with  the  Nemark  direct  integration  technique  [12].  The  closed  form 
solution  presented  herein  is  essentially  an  extension  of  the  static 
solution  presented  in  [11]  to  transient  solution.  The  finite  element 
model  and  the  closed-form  solution  are  based  on  a  shear-deformation 
theory  (see  Whitney  and  Pagano  [13])  of  layered  composite  plates,  which 
does  not  account  for  delamination  between  layers. 


Theory  and  Formulation 

The  displacement  field  in  all  of  the  simple  Timoshenko-type  shear- 
deformation  theories  are  based  on  the  following  displacement  field: 
Ujfx.y.z.t)  =  u(x,y,t)  +  z^fx.y.t) 


u2(x,y,z,t)  =  v(x,y,t)  +  z<t,y(x,y,t) 
u3(x,y,z,t)  =  w(x,y,t). 


0) 


Here  t  is  the  time,  u^  (i  =  1,2,3)  is  the  displacement  in  x-j -coordinate 

direction  (xj  =  x,  x2  =  y,  X3  =  z),  (u,v,w)  are  the  associated  midplane 

displacements,  and  <j>  and  4,  are  the  bending  slopes  in  the  xz  and  yz 

x  y 

planes.  The  coordinate  system  is  chosen  such  that  the  xy-plane 
coincides  with  the  midplane  of  the  plate. 

The  equations  of  motion  in  the  small -deflection  theory  of 
moderately  thick  plates  are  given  by 

Nl,x  +  N6,y  =  Pu,tt  +  R*x,tt 
N6,x  +  N2,y  ‘  Pv,tt  +  R<ly,tt 


°l,x  +  °2,y  *  Pw,tt  + 


(2) 
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Vx  +  M6,y  •  °1  =  l4x,tt  +  Ru,tt 
M6,x  +  M2,y  °2  l4y,tt  +  Rv,tt 

where  P,  R  and  I  are  the  normal,  coupled  normal -rotary,  and  rotary 
inertia  coefficients, 

h/2  9  zm+l  , 

(P.RJ)  =  /h/2  (l,z,zz)pdz  =  E  J  (l,z,zz)pWdz  (3) 
'  m  z 

m 

being  the  material  density  of  the  m-th  layer,  Ni,  Q-j,  and  Mi  are 
the  stress  and  moment  resultants. 

The  main  difference  between  the  usual  plate  theory  and  the 
laminated  plate  theory  is  reflected  in  the  plate  constitutive 
equations.  For  bimodular-material  plates  these  constitutive  equations 
are  given  by 


,  (1,j  =  1,2,6) 


e. 

J 

Dtj. 

Ki 

k2s 

k4^44 

k4k5S 

k4k5S45 

kZ$ 

*5*55 

Here,  Aij,  Bij,  D-j j ,  and  Sij  denote  the  respective  In-plane,  bending¬ 
stretching  coupling,  bending  or  twisting,  and  thickness-shear 
stiffnesses  defined  as  follows: 

h/2  9 

<W0u>  ■  L''-2-2  'W2  (,,j '  ',2,6) 


sii  ■  V1Jkl *  <I,J'4'5) 


where  h  is  the  total  thickness  of  the  plate,  denotes  the  plane- 

stress  reduced  stiffness  (1,j  refer  to  the  position  in  the  compliance 

matrix;  k  refers  to  the  sign  of  the  fiber-direction  strain:  k  =  1, 

tensile  and  k  =  2,  compressive;  and  z  refers  to  the  layer  number),  e.  and  K . 

J  J 

are  the  strains  and  curvatures  associated  with  the 


displacements  in  (1),  and  k^  are  the  shear  convection  coefficients. 


Closed-Form  Solutions 

In  [8,10,11]  it  was  shown  that  the  closed-form  solution  can  be 
derived  for  a  freely  supported,  lamianted,  bimodular,  rectangular  plate 
subjected  to  a  sinusoidally  distributed  normal  pressure.  Guided  by 
these  results,  we  consider  the  same  problem  but  subjected  to  sinusoidal 
distributed,  time-dependent  load, 

q(x,y,t)  =  (qQ  sin  ax  sin  py)T(t),  a  =  n/a,  p  =  u/b  (7) 
where  T(t)  is  a  known  function  of  time,  and  a  and  b  are  the  planform 
dimensions  of  the  plate.  The  boundary  condtions  are  given  by 
at  x  =  0,a:  v  =  w  =  =  Nj  =  =  0 

at  y  *  o,b:  u  =  w  =  4»x  =  =  0.  (8) 

For  the  loading  and  boundary  conditions  given  above,  the  governing 
equations  (2)  are  satisfied  exactly  (for  any  t  >  0)  by  the  following 
form  of  the  generalized  displacements: 

u(x,y,t)  =  U(t)^1(x,y) 
v(x,y,t)  =  V(t)<>2(x,y) 

w(x,y,t)  =  W(t)4>3(x,y)  (9) 

<l»x(x,y,t)  =  X(t)*1(x,y) 

4>y(x,y,t)  -  Y(t)4>2(x,y) 


where 
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♦j  =  cos  ax  sin  py,  ^  =  sin  cos  8Y>  =  sin  <*  sin  py.  (10) 

The  values  of  the  coefficients  U(t),  V(t),  W(t),  X(t)  and  Y(t)  are 
obtained  by  solving  the  following  ordinary  differential  equations  in 
time: 


£ 


[M] 


u 

u 

0 

V 

V 

0 

M 

+  [C] 

w 

= 

T 

X 

X 

0 

Y 

Y 

0 

(11) 


The  elements  of  [C]  and  [M]  are  given  in  [14]. 

It  should  be  pointed  out  that  the  plate  stiffnesses  A^j,  B^j,  and 

D^j  depend  on  the  neutral -surface  positions  Znx  and  Zny  associated  with 

the  fiber-direction  strains  e  and  e  ,  respectively.  The  neutral- 

x  y 

surface  locations  are  given  by. 


znx(t)  -  -U(t)/X(t) 

Zny(t)  =  -V(t)/Y(t). 


(12) 


which  are  constants  for  any  fixed  time.  Thus,  for  any  t  >  0,  the 
neutral  surfaces  are  planes,  as  In  the  case  of  static  bending.  For  the 
detailed  computations  of  A^,  B^,  and  D^.,  the  reader  Is  referred  to 
[8,9,11]. 


Next,  we  consider  the  numerical  solution  of  (11)  by  the  Newmark 
direct  Integration  technique.  Equation  (11)  is  of  the  general  form 

[M]{a}  +  [K]  {a}  =  (F)  (13) 

where  (a)  denotes  the  column  of  the  generalized  displacements.  In 
Newmark's  method,  the  solution  and  its  derivative  with  respect  to  time 
are  approximated  by 


I 
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{A}n+1  =  {A}n  +  At{A}n  +  i(\  -  p){A}n  +  p{A)n+1](At)2 

{^}n+1  =  +  ^0  ~  a^A^n  +  a^AWl^At  04) 

where  a  and  p  are  the  parameters  that  control  the  stability  of  the 
scheme,  {A}n  denotes  the  value  of  { A }  at  time  t  =  n  At.  The  choice 
a  =  0.5  and  p  =  0.25  is  known  to  give  an  unconditionally  stable  scheme 
for  linear  problems. 

Using  of  the  approximations  (14)  in  (11)  and  rearranging  the  terms, 
we  arrive  at 

KJW„+1  ■  (f  05) 

where 

[K]  =  [K]  +  aQ[M] 

(Ff  =  {Ffn+1  +  ^M^ao^A^n  +  al^A^n  +  a2^AU  06) 

aQ  =  l/(pAt)  ,  ax  =  aQAt  ,  a2=-£jj-l. 

Starting  with  initial  values  of  { a}0 ,  {a}0  and  {a}0,  equation  (15)  can 
be  solved  repeatedly  for  (a)  at  successive  values  of  time.  The  values 
of  {A}n+1  and  { A}n+i  can  be  computed  from  (14). 

Finite-Element  Formulation 

As  pointed  out  earlier,  the  finite  element  model  used  in  the 
present  study  is  the  same  as  that  employed  in  [9].  We  shall  not  repeat 
the  formulation  here,  but  only  point  out  the  additional  steps  involved 
in  the  transient  analysis.  The  finite-element  model  in  the  present  case 
results  in  the  following  equation  for  an  element 

[M]{'M  +  CK]{a}  =  {F},  (17) 

where  {a}  denotes  the  column  of  the  nodal  values  of  the  generalized 
displacements.  The  elements  of  the  mass  matrix  [M]  and  the  stiffness 


matrix  [k]  are  given  in  [14,15].  Equation  (17)  is  integrated  using  the 
direct  integration  technique  described  above. 

It  should  be  noted  that  no  restriction  is  placed  on  the  loading  and 

boundary  conditions  In  the  finite-element  analysis.  When  a  uniform 

loading  is  used,  the  neutral -surface  locations  are  not  independent  of 

the  location  (x,y),  and  therefore,  the  expressions  for  Ay,  By  and  Dy 

derived  for  the  closed-form  solutions  are  not  valid  in  the  entire  plate 

in  the  finite-element  analysis.  Since  the  coefficients  K. .  are 

*  J 

evaluated  at  the  Gauss  points,  the  stiffnesses  Ay,  By  and  Dy  are  also 
evaluated  at  the  Gauss  points  by  using  the  neutral -surface  positions 
computed  at  the  Gauss  points.  This  is  equivalent  to  the  assumption  that 
the  neutral  surfaces  are  element-wise  bilinear. 

In  the  actual  calculation  of  the  stiffnesses,  whether  for  the 
closed-form  solution  or  In  the  finite-element  analysis,  the  neutral- 
surface  locations  are  not  known  a-priori .  Therefore  an  iterative 
technique  must  be  employed  to  compute  the  neutral -surface  locations. 

The  iteration  procedure  begins  with  assumed  values  of  Znx  and  Zny  (say, 
Znx  =  Zny  =  0)  and  then  the  stiffnesses  Ay  etc.  are  computed  using 
these  values.  After  obtaining  the  generalized  displacements,  the 
neutral -surface  locations  are  recomputed.  Using  these  new  values  of  Znx 
and  Zny,  the  stiffnesses  for  the  next  iteration  are  computed.  This 
procedure  is  repeated  until  the  difference  between  any  two  consecutive 
values  of  Znx  (and  Zny)  differ  by  a  small  preselected  value  (say  0.1%). 

Numerical  Results  and  Discussion 

In  the  following,  numerical  results  are  presented  for  rectangular 
plates  made  of  two  bimodular  materials:  aramld  cord-rubber  (AR)  and 
polyester  cord-rubber  (PR),  which  are  used  in  automobile  tires.  The 
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material  properties  for  these  two  materials  are  given  in  [7-11,14],  and 

are  listed  in  Table  1  again  for  convenience.  In  the  finite  element 

method,  a  2  x  2  mesh  of  nine-node  Isoparametric  elements  in  the  quarter 

2 

plate  was  used.  The  shear  correction  coefficients  kj  were  chosen  to  be 
5/6. 


Table  1  Material  properties  for  aramid  cord-rubber  and  polyester  cord- 
rubber,  unidirectional,  bimodulus  composite  materials. 


Property 

Aramid 

-Rubber 

Polyester  Rubber 

Tensile 

Compressive 

Tensile 

Compressive 

En  (GPa) 

3.58 

0.012 

0.617 

0.0369 

e22  (Gpa) 

0.00909 

0.012 

0.008 

0.0106 

v,, 

0.416 

0.205 

0.475 

0.185 

(l]2  *  Git  (GPa) 

0.0037 

0.0037 

0.00262 

0.00267 

G„  (GPa) 

0.0029 

0.00499 

0.00233 

0.00475 

P  (Gpa) 

1.0 

1.0 

1.0 

1.0 

Due  to  the  lack  of  other  results  in  the  literature,  comparisons  are 
made  between  the  present  closed-form  solution  (CFS)  and  finite  element 
solutions  (FES).  The  selection  of  the  time  step  was  guided  by  the 
stability  criterion  given  by  Tsui  and  Tong  [16]  for  moderately  thick 
plates  (see  [15]). 

First,  the  effect  of  the  time  step  on  the  transverse  deflection  and 
the  neutral -surface  location  of  a  single-layer  (0°)  square  plate  (a  =  b 
=  1,  a/h  =  5,  aramid-rubber)  under  suddenly  applied  step  loading  (with 
q0  =  1  in  eqn.  (7))  was  studied.  Table  2  contains  the  results  for  the 
transverse  deflection  and  the  neutral -surface  location  for  various  time 
steps.  Note  that  the  solutions  obtained  using  At  =  0.5  differ  from 
those  obtained  using  At  =  0.05  by  only  2%. 

Next,  the  finite-element  results  are  validated  by  comparing  with 
the  closed-form  solutions  of  a  two-layer  (0°/90°)  square  plate  (a  =  b  = 
1,  a/h  =  5,  aramid-rubber)  under  suddenly  applied  step  loading  (7). 


I 


Table  2  Effect  of  the  time  step  on  the  transverse  deflection  and  the  neutral -surface  location 
of  a  single-layer  (aramid-rubber)  square  plate  under  suddenly  applied,  sinusoidally 
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Table  3  contains  the  results  for  the  transverse  deflection,  transverse 
velocity,  and  the  neutral  surface  location,  znx.  The  finite  element 
results  are  in  excellent  agreement  (less  than  one  percent  error)  with 
the  closed-form  solution.  Both  the  closed-form  and  finite-element 
solutions  were  obtained  using  At  =  0.5,  and  q0  =  1  in  eqn.  (7).  In  view 
of  the  close  agreement  between  the  closed-form  and  finite-element 
solutions,  and  in  the  interest  of  saving  the  computational  time  (the 
finite  element  method  took  30  min.  of  CPU  time  compared  to  5  sec.  for 
the  closed  form  solution  presented  in  Table  3),  most  of  the  results  to 
be  discussed  were  obtained  using  the  closed-form  method.  Of  course,  for 
the  uniformly  distributed  load  case,  the  finite  element  method  was 
empl oyed. 

Figure  1  contains  the  plots  of  the  transverse  deflection  and  the 
neutral -surface  location  versus  time  for  a  single-layer  (0°),  square 
plate  (a  ■  b  =  1,  a/h  =  5,  aramid-rubber)  subjected  to  suddenly  applied 
step  loading  (q0  =1,  At  =  0.5).  The  figure  shows  the  effect  of  the 
shear  deformation  and  the  bimodular  action  on  the  amplitude  and  period 
of  the  solutions.  Clearly,  the  effect  of  both  the  shear  deformation  and 
the  bimodular  action  is  very  pronounced  and  hence  cannot  be  neglected  in 
the  analysis.  It  should  be  pointed  out  that  the  value  of  Znx/h  for  the 
time  interval  8  <  t  <  36,  coincides  with  the  Znx/h  of  the  associated 
static  case  (see  [11]). 

To  further  Investigate  the  effect  of  the  plate  thickness  on  the 
dynamic  response,  the  same  problem  (as  In  Figure  1)  was  solved  using  a/h 
=10.  The  nondimensional Ized  transverse  deflection,  w^wh^^/q  a£, 
versus  time  Is  shown  in  Fig.  2.  The  effect  of  the  thickness  on  the 
amplitude  and  period  of  the  deflection  is  apparent  from  the  figure. 


11 


Table  3  Comparison  of  the  neutral -surface  location,  transverse 
deflection  and  transverse  velocity  as  obtained  by  the 
closed-form  solution  (CFS)  and  the  finite  element  method 
(FEM)  for  a  two-layer  (0°/90°)  square  plate  (aramid- 
rubber)  under  suddenly  applied  transverse  step  loading'. 


t 

— _ ?nx 

/h 

W  X 

10-3 

W  X 

10'2 

CFS 

FEM 

CFS 

FEM 

CFS 

FEM 

4 

0.3780 

0.3804 

0.0330 

0.0330 

0.1674 

0.1671 

8 

0.4387 

0.4385 

0.1277 

0.1277 

0.2959 

0.2966 

12 

0.4377 

0.4370 

0.2597 

0.2596 

0.3486 

0.3477 

16 

0.4395 

0.4396 

0.3904 

0.3900 

0.2908 

0.2902 

20 

0.4410 

0.4402 

0.4832 

0.4826 

0.1614 

0.1603 

24 

0.4388 

0.4389 

0.5127 

0.5112 

-.0208 

-.0223 

28 

0.4419 

0.4411 

0.4684 

0.4665 

-.1912 

-.1920 

32 

0.4373 

0.4372 

0.3656 

0.3632 

-.3116 

-.31 36 

36 

0.4412 

0.4410 

0.2307 

0.2284 

-.3457 

-.3438 

.40 

0.4279 

0.4255 

0.1044 

0.1026 

-.2745 

-.2735 

a mJteb 
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Figure  1.  Effect  of  the  transverse  shear  deformation  and  bimodulus  action 
on  the  transient  response  of  single-layer  (0°)  square  plates 
(a  3  b  =  1,  h  =  0.2,  q  =  1 ,  At  =  0.5)  of  aramid-rubber  under 
suddenly  applied,  sinusoidally  distributed,  step  loading. 
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The  transverse  deflection,  transverse  velocity,  and  the  neutral  - 
surface  location  versus  time  for  a  single-layer  (0°),  square  plate  of 
polyester-rubber  material  (a/h  =  5,  q0  =  1,  At  =  0.5)  under  suddenly 
applied  step  loading  are  shown  in  Fig.  3.  The  figure  also  contains  the 
plot  of  the  transverse  deflection  versus  time  for  a  rectangular  plate 
(a/b  =  2;  everything  else  is  the  same  as  in  the  square  plate). 

In  Fig.  4,  results  of  two-layer  (0°/90°)  square  plates  (a  =  b  =  1, 
a/h  =  5,  At  =  0.5,  q0  =  1)  of  aramid-rubber  and  polyester-rubber  are 
presented.  Once  again,  the  excellent  agreement  between  the  closed-form 
solution  and  the  finite-element  solution  is  observed.  Note  also  that 
the  neutral -surface  locations,  both  Znx/h  and  Zny/h,  for  the  polyester- 
rubber  plate  are  bounded  above  by  Znx/h  and  below  by  Zny/h  of  the 
aramid-rubber  plate. 

Figure  5  contains  plots  of  solutions  of  a  single-layer  (0°)  and  a 
two-layer  (0°/90°)  square  plate  (a  =  b  =  1 ,  a/h  =  5,  At  =  0.5,  q0  =  1 ) 
of  aramid-rubber  under  suddenly  applied  impulse  loading, 
q(x,y,t)  =  q0sin(^-)sin(^)H(t  -  5). 
where  H(t)  is  the  Heavyside  step  function.  Note  that  the  neutral - 
surface  location,  Znx/h,  does  not  reach  a  steady  value  like  In  the  case 
of  step  loading,  but  the  mean  value  is  about  the  same  as  that  in  the 
case  of  step  loading. 

Thus  far,  the  transverse  load  was  assumed  to  be  sinusoidal  (with 
respect  to  x  and  y)  to  facilitate  closed-form  solutions.  Since  the 
finite  element  method  does  not  have  any  limitation  on  the  load 
distribution,  one  can  use  the  finite  element  model  developed  herein  for 
the  analysis  of  plates  subjected  to  uniformly  distributed  loads.  Figure 
6  contains  the  plots  of  the  transverse  deflection,  transverse  velocity 


Figure  3.  Transient  response  of  a  single-layer  (0°),  bimodular,  rectangular 
plate  of  polyester-rubber  material  under  suddenly  applied,  sinu¬ 
soidally  distributed,  step  loading. 


o  Aramid-rubber-FES 

—  Aramid-rubber-CFS 

—  Polyester-rubber-CFS 


w  x  10 
(Zn/h)  0.1 


V  / 

•  A  ' 


4  8  12  16  20  24  29  32  36  40  44  48  52  t 


Figure  4.  Transient  response  of  two-layer  (0°/90°),  bimodular,  square 
plates  (a  *  b  *  1,  h  =  0.2,  q0  =  1 ,  at  =  0.5)  of  aramid- 
rubber  and  polyester-rubber  materials  subjected  to  suddenly 
applied,  sinusoidally  distributed,  step  loading. 


(Z  /h)/40 


and  the  neutral -surface  location,  Znx/ h,  for  a  single-layer  (0°)  and  a 
two-layer  (0°/90°)  square  plate  (a  =  b  =  1 ,  a/h  =  5,  At  =  0.5,  q0  =  1 ) 
of  aramid-rubber  under  suddenly  applied  step  loading. 


q(x,y,t)  =  qQH(t). 

The  response  curves  resemble  those  of  ordinary  plates  (see  [15]). 
Although  the  amplitude  of  the  transverse  deflection  is  quite  large 
compared  to  that  of  a  sinusoidally-loaded  plate,  the  magnitude  of  the 
neutral -surface  location  is  about  the  same.  This  completes  the 
discussion  of  the  numerical  examples. 

Summary  and  Conclusions 

The  transient  analysis  of  bimodular,  composite,  rectangular  plates 
is  presented.  Finite  element  as  well  as  closed-form  solutions  are 
presented  for  rectangular  bimodular  plates  of  aramid-rubber  and 
polyester-rubber  materials.  The  finite  element  solutions  are  found  to 
agree  very  well  with  the  closed-form  solutions.  From  the  dynamic 
(transient)  response  of  the  bimodular  plates,  it  is  apparent  that  the 
shear  deformation  and  plate  thickness  (irrespective  of  the  transversse 
shear  strains)  increase  the  amplitude  and  period  of  the  transverse 
deflection.  The  results  of  the  present  study  should  be  of  interest  to 
composite-materials  designers  and  researchers. 
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